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ABSTRACT 

This  paper  examines  the  multi-level,  primitive  equation,  global  ocean  circulation 

model  of  Semtner  and  Chervin  for  its  ability  to  simulate  the  seasonal  cycle  in  the  tropical 
Pacific  Ocean.  The  result  of  a  20-year  integration  of  this  model  using  annual  mean  wind 
forcing  was  reported  in  Semtner  and  Chervin  (1988).  This  was  the  first  global  eddy- 
resolving  ocean  calculation  and  it  showed  many  realistic  features  of  ocean  circulation. 
The  phase  of  the  simulation  analyzed  in  this  report  incorporates  seasonally  varying  wind 
forcing  from  the  Hellerman  and  Rosenstein  (1983)  global  data  set.  These  wind  stress 
values  were  defined  on  a  grid  with  2°  spacing  which  have  been  interpolated  to  the  one- 
half  degree  grid  points  of  the  Semtner  and  Chervin  model.  There  is  no  interannual 
variability  in  the  wind  fields  of  this  data  set.  The  results  presented  here  are  from  the 
fourth  year  of  a  10-year  seasonal  cycle  run. 

The  upper  oceanic  circulation  of  the  equatorial  Pacific  Ocean  consists  of 
alternating  westward  and  eastward  flows  that  vary  seasonally  in  strength  and  extent.  The 
eastward  flows  investigated  in  this  study  are  the  Equatorial  Undercurrent  and  North 
Equatorial  Countercurrent;  the  westward  currents  are  the  North  and  South  Equatorial 
Currents.  The  Equatorial  Undercurrent  is  most  intense  late  in  the  Northern  Hemisphere 
spring  and  is  weak  in  the  late  fall  at  which  time  the  North  Equatorial  Countercurrent  and 
the  South  Equatorial  Current  are  most  intense.  This  result  is  in  agreement  with 
observations  and  other  model  studies. 
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I.  INTRODUCTION 

A.   BACKGROUND 

Ocean  numerical  modelers  for  many  years  have  been  limited  in  the  scope  of  the 
ocean  model  studies  they  can  achieve  because  of  limitations  in  computer  power  and 
resources.  Oceanographers  have  watched  numerical  modeling  in  the  meteorology 
community  grow  into  a  critical  tool  for  understanding  and  predicting  the  weather.  In 
contrast  to  atmospheric  modeling,  world  ocean  "climate"  modeling  requires  a  thousand- 
fold increase  in  computer  power  due  to  the  order  of  magnitude  increase  in  each 
horizontal  direction  required  to  model  oceanic  features  versus  atmospheric  (because  the 
Rossby  radius  of  deformation  at  midlatitudes  is  typically  400  km  in  the  atmosphere  and 
40  km  in  the  ocean).  Integration  time  for  numerical  ocean  models  requires  another  ten- 
fold increase  due  to  the  longer  time  scales  of  oceanic  phenomenon.  Oceanographers  also 
suffer  from  a  paucity  of  observational  data  with  which  to  initialize  models.  Because  of 
these  factors,  oceanographers  have  lagged  behind  their  colleagues  in  meteorology  in 
many  phases  of  applications,  research  with  and  prognostic  use  of  numerical  modeling. 

Modern  generations  of  supercomputers  are  finally  available  to  oceanographers  that 
are  powerful  enough  to  run  "full  physics"  general  circulation  simulations  on  a  global 
scale.  Still  the  resolution  of  these  ocean  general  circulation  models  are  marginal  in  their 
ability  to  resolve  mesoscale  eddies,  the  oceanic  equivalent  of  atmospheric  cyclones  and 
anticyclones.    A  computer  simulation  of  the  tropical  oceans  holds  the  most  promise  in 
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having  "realistic"  currents  and  thermal  structure  from  a  marginal  eddy-resolving  global 
ocean  model  because  of  the  larger  Rossby  radius  of  deformation  (-140  km)  near  the 
equator.  This  paper  is  an  analysis  of  the  results  in  the  tropical  Pacific  Ocean  of  one  such 
global  simulation. 

In  recent  years,  a  large  amount  of  public  attention  has  been  focused  on  the  question 
of  "global  warming."  The  growing  global  environmental  movement  owes  its  popularity 
in  large  part  to  public  perceptions  and  concerns  for  this  issue.  Vast  sums  of  research 
dollars  are  being  funneled  into  investigating  climate  change.  Scientists  have  long  known 
that  the  earth  undergoes  natural  climate  change  over  long  periods.  The  question  is 
whether  or  not  human  activities  since  the  industrial  revolution  have  accelerated  climate 
change,  change  which  may  have  unforseen  and  possibly  catastrophic  consequences  for 
much  of  the  world's  population  in  the  form  of  starvation,  social  and  political  upheaval. 

The  ocean  is  known  to  be  a  sink  for  carbon  dioxide  and  can  store  and  redistribute 
a  vast  amount  of  heat.  The  tropical  oceans  play  a  critical  but  poorly  understood  role  in 
this  redistribution  of  heat.  The  feedback  between  the  atmosphere  and  the  ocean  will 
require  the  most  sophisticated  coupled  ocean-atmosphere  models  in  order  to  accurately 
simulate  climate  variability. 

Scientists  cannot  hope  to  have  an  accurate  coupled  ocean-atmosphere  model  for 
studies  of  long-term  climate  change  until  we  first  have  an  understanding  and  can  model 
oceanic  heat  flux  and  global  ocean  circulation.  The  Semtner  and  Chervin  global 
circulation  model,  hereafter  referred  to  as  the  SC-GCM,  is  an  attempt  to  model  the 
global  ocean  circulation. 


The  SC-GCM,  is  a  one-half  degree  primitive  equation  model  that  ran  for  a  20-year 
integration  using  annual  mean  wind  forcing.  This  simulation  showed  many  realistic 
features  of  oceanic  circulation  (Semtner  and  Chervin,  1988).  The  next  step  in  the 
evolution  of  this  model  was  to  incorporate  seasonally  varying  forcing  of  wind  from  the 
Hellerman  and  Rosenstein  (1983)  global  data  set  and  seasonal  heat  and  salt  fluxes  from 
the  Levitus  (1982)  data  base.  Assuming  the  forcing  is  reasonably  correct,  the  simulation 
should  produce  a  strong  seasonal  cycle  in  the  tropical  oceans.  The  results  of  the  model 
presented  are  from  the  fourth  year  of  a  10-year  seasonal  cycle  run. 

B.      PACIFIC  OCEAN 

The  upper  oceanic  circulation  of  the  equatorial  Pacific  Ocean  consists  of  alternating 
westward  and  eastward  flows  as  shown  schematically  in  Figure  1.1.  The  four  major 
currents  consist  of  three  surface  currents,  the  westward  flowing  North  Equatorial  Current 
(NEC)  extending  between  latitudes  10°N  and  20°N,  the  westward  flowing  South 
Equatorial  Current  (SEC)  between  3°N  and  10°S,  and  between  these  two  westward 
surface  currents,  the  eastward  flowing  North  Equatorial  Countercurrent  (NECC).  The 
Equatorial  Undercurrent  (EUC)  is  the  fourth  major  current  and  is  a  submerged  current 
flowing  eastward  along  the  equator  from  2°N  to  2°S.  Observationalists  have  documented 
that  the  Equatorial  Undercurrent  is  most  intense  during  the  Northern  Hemisphere  spring 
and  is  weak  in  the  fall  at  which  time  the  North  Equatorial  Countercurrent  and  the  South 
Equatorial  Current  are  most  intense.  The  existence  of  these  four  currents  has  been 
known  for  many  years,  but  scientists  are  still  far  from  a  comprehensive  understanding 


of  their  variability  both  seasonally  and  interannually.  Oceanographers  know  in  a 
qualitative  sense  that  these  currents  vary  seasonally  and  spatially  in  strength  and  vertical 
and  horizontal  extent,  i.e. ,  volume  transport.  They  have  no  clear  understanding  of  their 
climatological  state  on  an  annual  and  seasonal  basis.  Because  of  the  large  interannual 
variability  in  the  forcing  over  the  tropical  Pacific  Ocean,  such  a  climatology  would 
require  observations  over  many  years  and  over  a  huge  expanse  of  ocean.  General 
circulation  numerical  ocean  models  are  important  tools  in  assisting  oceanographers  to 
arrive  at  this  understanding. 

The  organization  of  this  thesis  is  as  follows.  A  brief  description  of  the  model  used 
occurs  in  Chapter  II ,  along  with  the  trade  wind  system  which  drives  the  surface  currents 
in  the  tropics.  Next,  the  major  equatorial  currents  in  the  model  are  examined  in  the 
model  and  compared  to  observational  studies.  Additionally,  the  seasonal  variations  of 
these  currents  and  the  associated  variation  of  wind  forcing  are  examined  (Chapter  III). 
Several  lesser  known  off-equatorial  currents  are  also  resolved  in  the  SC-GCM  and  are 
discussed.   Finally,  conclusions  and  recommendations  are  presented  in  Chapter  IV. 


H.  GLOBAL  OCEAN  MODEL 

The  model  used  in  this  study  is  the  multi-level,  primitive  equation,  robust- 
diagnostic  global  ocean  circulation  model  of  Semtner  and  Chervin  (1988).  This  model 
has  been  described  in  detail  in  the  existing  oceanographic  literature.  What  follows  then 
is  a  general  description  of  the  global  ocean  model  with  the  exception  that  the  monthly 
varying  wind  stress  field  in  the  Pacific  basin  used  in  the  model  calculation  will  be 
discussed  thoroughly.  For  the  interested  reader  additional  information  on  the  SC-GCM 
formulation  can  found  in  Semtner  (1986  a,b). 

A.      MODEL  PARAMETERS 

The  domain  of  the  SC-GCM  is  all  longitudes  from  latitudes  75 °S  to  65 °N.  The 
grid  spacing  is  1/2°  degree  in  the  horizontal  and  there  are  20  levels  in  the  vertical. 
Realistic  coasts  and  bathymetry  (slightly  smoothed)  are  used. 

Wind  stress  is  specified  from  Hellerman  and  Rosenstein  (1983)'s  monthly  averaged 
climatological  analysis  of  ship  winds.  This  global  data  set  defines  values  on  a  grid  of  2° 
spacing.  For  use  in  the  1/2°  spacing  of  the  SC-GCM,  the  wind  data  have  been 
interpolated  to  the  grid  points  of  the  model.  A  smooth  transition  curve  was  applied  to 
the  mean  monthly  values  of  wind  stress  from  one  month  to  the  next  to  avoid  "jolting" 
the  model  with  sudden  and  unrealistic  climatological  changes.  It  should  be  emphasized 
that  at  this  stage  of  the  simulation,  the  wind  stress  used  has  seasonal  but  no  interannual 


variability.  Thus  the  SC-GCM  is  not  designed  to  model  interannual  phenomenon  in  the 
Pacific  Ocean  such  as  the  El  Nino/Southern  Oscillation  (ENSO). 

Vertical  mixing  in  the  model  uses  the  Richardson-number  dependent 
parameterization  of  Pacanowski  and  Philander  (1981).  Effects  of  eddies  smaller  than  the 
eddies  resolved  by  the  model  are  parameterized  by  biharmonic  mixing.  The  SC-GCM 
is  initialized  with  temperature  and  salinity  values  from  the  Levitus  (1982)  global  data 
base.  The  upper  25  meters  of  the  water  column  in  the  model  are  forced  to  these  Levitus 
values  on  a  monthly  time  scale  in  the  method  of  Haney  (1971).  This  serves  as  a  proxy 
for  heat  and  moisture  fluxes.  Below  700  meters  temperature  and  salinity  values  are 
restored  to  Levitus  data  on  a  three-year  time  scale.  Restoring  within  the  thermocline  (25- 
700  m)  occurs  only  in  latitude  bands  65°S-75°S  and  55°N-65°N,  also  on  a  three-year 
time  scale.  In  all  other  latitude  bands,  including  the  Pacific  equatorial  region  of  this 
study,  the  thermocline  is  unconstrained  and  free  to  resolve  its  own  temperature  and 
salinity  structure,  currents  and  eddies.  The  model  runs  with  a  900  second  time  step 
between  iterations. 

The  model  geometry  has  been  simplified  somewhat  by  connecting  most  major 
islands  to  continental  land  masses  and  shoaling  remaining  islands  to  a  depth  of  100 
meters.  For  the  Pacific  portion  of  the  global  model  (Figure  2.1),  New  Guinea  is 
connected  to  Australia  and  the  western  Pacific  islands  of  Celebes,  Mindanao  and 
Halmahera  are  submerged.  Borneo  and  much  of  the  Philippines  are  joined  to  the  East 
Asian  land  mass  and  the  marginal  seas  in  between  are  filled  in.  All  other  islands  in  the 
Pacific  are  submerged  with  the  exception  of  New  Zealand.    As  will  be  shown  in  the 


analysis  of  the  results,  this  simplified  geometry  and  bathymetry  in  the  western  Pacific  has 
a  significant  impact  on  the  model  results  for  this  area.  Bathymetry  in  the  model  is 
determined  by  the  number  of  vertical  grid  boxes  at  each  location  up  to  the  limit  of  20 
(maximum  model  depth  of  5200  meters).  The  boxes  are  stacked  from  the  surface  down 
until  the  local  bottom  is  reached  at  each  horizontal  grid  point.  The  20  vertical  grid  levels 
have  varying  thicknesses  to  allow  concentration  of  boxes  in  the  upper  portion  of  the 
ocean  for  enhanced  thermocline  resolution.  The  bathymetry  of  the  model  in  the  Pacific 
is  shown  in  detail  and  at  various  depths  in  Figure  2.2. 

The  SC-GCM  has  completed  32.5  years  of  simulation  with  the  first  22.5  years 
using  annual  mean  wind  forcing,  heat  and  salt  fluxes.  Results  from  the  annual  mean  run 
were  discussed  in  Semtner  and  Chervin  (1988).  Monthly  mean  forcing  has  been  run  for 
the  remaining  10  years.  At  the  end  of  each  three-day  period  of  simulation,  a  "snapshot" 
was  taken  of  the  global  ocean.  Years  20-30  of  the  simulation  are  archived  at  these  three- 
day  intervals  for  ongoing  analysis. 

B.      MONTHLY  MEAN  WIND  STRESS 

The  equatorial  currents  in  the  Pacific  are  driven  by  the  trade  wind  system  that 
persists  throughout  the  year.  The  variations  in  the  wind  stress  pattern  cause  variations 
in  the  position  of  the  Intertropical  Convergence  Zone  (TTCZ),  the  area  between  the 
northern  and  southern  hemispheric  trade  wind  belts  where  their  convergence  occurs.  The 
seasonal  variation  in  the  wind  stress  field  does  not  rival  the  monsoonal  reversal  seen  in 
the  Indian  Ocean,  but  some  monsoonal  effects  are  present  in  the  western  Pacific  Ocean. 


Because  the  surface  currents  are  wind-driven,  an  understanding  of  the  major  currents  of 
the  equatorial  Pacific  requires  an  understanding  of  the  equatorial  wind  field.  Figures 
2.3-2. 14  are  vector  diagrams  of  the  annual  cycle  of  wind  stress  in  the  equatorial  Pacific 
Ocean  from  the  Hellerman  and  Rosenstein  (1983)  global  data  set.  The  wind  stress 
pattern  will  be  described  for  the  western  and  central  equatorial  Pacific  first,  then  for  the 
simpler  pattern  of  the  eastern  Pacific. 

In  January  (Figure  2.3a),  strong  northeasterly  wind  flow  is  prevalent  in  the  central 
and  western  Pacific  from  the  atmospheric  pressure  gradient  around  the  Asiatic  high.  The 
northeasterly  flow  extends  to  the  equator  in  the  central  Pacific.  In  the  western  Pacific 
the  wind  field  is  more  complicated  as  it  becomes  northerly  and  even  northwesterly  as  it 
turns  along  the  Philippines  and  New  Guinea  coast.  The  southeasterly  (SE)  trades  are 
weak  and  are  evident  along  the  18-20°S  latitudinal  band.  The  ITCZ  is  the  broad  and 
diffuse  area  of  weak  winds  along  10°S  in  the  western  Pacific  and  extending  from  the 
equator  to  approximately  16°S  in  the  central  Pacific.  The  northeasterly  (NE)  trades  are 
at  their  strongest  during  this  month  with  maximum  wind  stress  values  of  approximately 
1.9  dyn  cm"2  in  the  Philippine  Sea.  The  southeast  trades  are  much  weaker,  with 
maximum  values  of  0.9  dyn  cm"2  near  the  northeastern  coast  of  Australia.  The  SE 
trades  weaken  toward  the  east. 

The  February  wind  pattern  (Figure  2.3b)  is  little  changed  from  January.  The 
northeast  monsoon  and  its  associated  winds  have  weakened  slightly. 

The  azimuthal  pattern  of  wind  stress  in  March  (Figure  2.4a)  is  nearly  identical  to 
the  first  two  months  except  that  the  northeast  monsoon  has  weakened  with  wind  stress 
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values  of  1.2  dyn  cm"2  in  the  western  equatorial  Pacific.  The  southeast  trades  have 
increased  in  strength,  but  they  remain  as  a  broad  area  of  weak  winds  extending  from  the 
equator  southward. 

In  April  (Figure  2.4b),  there  is  a  pronounced  weakening  of  the  NE  trades  beginning 
in  the  western  Pacific  Ocean  as  the  Asiatic  high  breaks  down.  Wind  stress  values  are 
a  maximum  of  0.7-0.8  dyn  cm"2.  In  the  central  Pacific  the  NE  trades  have  weakened 
only  slightly.  The  SE  trade  winds  have  now  increased  to  the  point  that  they  have  caused 
a  reversal  of  the  previous  southwesterly  flow  over  the  Indonesian  archipelago.  The  flow 
is  now  an  onshore  southeasterly  flow.  Additionally  the  SE  trades  have  begun  to  move 
the  ITCZ  north  from  5°S  to  10°S. 

May  shows  a  continuation  of  the  pattern  of  the  previous  month  (Figure  2.5a).  The 
NE  trades  continue  to  weaken.  They  are  more  easterly  in  the  central  Pacific  and  veer 
along  the  Philippine  coast  to  become  indistinguishable  from  the  SE  trades.  The  SE  trades 
have  not  yet  become  as  strong  as  the  NE  trades  in  the  central  Pacific  but  they  have 
changed  the  position  of  the  ITCZ  two-three  degrees  in  latitude  to  a  latitude  band  between 
2°S-7°S. 

In  June  (Figure  2.5b),  the  continuing  strengthening  of  the  SE  trades  in  the  western 
Pacific  Ocean  has  pushed  the  ITCZ  to  the  north  of  the  equator.  Cross-equatorial  flow 
from  south  to  north  is  now  evident  throughout  the  entire  central  and  western  tropical 
Pacific  Ocean  expanse.  The  strength  of  the  wind  stress  field  has  changed  only  slightly 
from  the  previous  month,  the  major  change  being  this  more  pronounced  cross-equatorial 
flow. 


By  July  (Figure  2.6a),  the  influence  of  the  summer  heat  low  over  Asia  has  resulted 
in  a  reversal  of  the  wind  flow  along  the  Philippine  coast  from  that  of  the  January  wind 
regime.  The  winds  in  the  western  tropical  Pacific  Ocean  now  have  a  southerly 
component  as  the  SE  trades  continue  to  strengthen  to  a  value  of  1 .5  dyn  cm-2.  Northeast 
of  Australia  the  influence  of  the  southeasterlies  has  pushed  the  ITCZ  to  5°N.  The  NE 
trades  in  the  central  Pacific  have  continued  to  weaken  and  are  approximately  1.0-1.2  dyn 
cm-2. 

The  pattern  in  August  (Figure  2.6b)  is  little  changed  from  the  previous  month.  The 
ITCZ  now  spans  latitudes  7°-9°N  in  both  the  central  and  western  Pacific  Ocean  basin. 

The  September  wind  stress  pattern  (Figure  2.7a)  continues  the  summer  pattern 
described  in  August.  There  is  however  a  noticeable  weakening  (.2-. 4  dyn  cm"2)  of  the 
SE  trades  in  the  western  Pacific.  The  NE  trades  have  also  weakened  (to  values  slightly 
above  1.0  dyn  cm-2)  and  the  ITCZ  is  now  at  10°N,  its  furthest  northward  extension  for 
this  region  of  the  Pacific. 

As  the  SE  monsoon  dies  out  in  October  over  east  Asia,  the  NE  trades  begin  to 
establish  themselves  again.  This  is  evident  over  the  western  Pacific  (Figure  2.7b)  as  the 
winds  begin  to  reverse  again  along  the  Philippine  coast.  There  is  still,  however,  cross- 
equatorial  flow  from  the  SE  trades.  The  ITCZ  is  nearer  the  equator  in  the  western 
Pacific  but  maintains  near  10°N  in  the  central  Pacific. 

In  November  (Figure  2.8a),  the  NE  trades  show  a  pronounced  increase  in  strength 
(0.5  dyn  cm'2)  over  the  Philippine  Sea.  The  SE  trades  have  begun  to  retreat  and  there 
is  now  a  broad  area  of  very  light  winds  from  7°S  to  7°N.    In  the  central  Pacific,  the 
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increase  in  the  NE  trades  is  also  pronounced  as  the  Asiatic  cold  high  begins  to  build. 
There  are  two  areas  of  convergent  winds  in  the  central  Pacific,  one  along  7°N  and  a 
second  but  weaker  one  along  10°S. 

In  December  (Figure  2.8b),  the  NE  trades  have  again  become  dominant  throughout 
the  western  and  central  equatorial  Pacific  Ocean.  The  NE  winds  turn  along  the  northern 
New  Guinea  coast  to  become  northwesterly.  The  pattern  of  two  areas  of  trade  wind 
convergence  in  the  central  Pacific  makes  the  identification  of  the  ITCZ  difficult,  there 
being  convergence  of  SE  and  NE  trades  along  both  5°N  and  12°S. 

Figures  2.9-2.14  illustrate  the  wind  stress  pattern  for  the  eastern  Pacific.  The 
annual  cycle  is  similar  to  the  rest  of  the  Pacific  but  the  pattern  is  simpler  and  more 
regular.  The  NE  trades  are  strongest  in  January  (Figure  2.9a)  and  February  (Figure 
2.9b)  when  the  SE  trades  are  weakest  as  in  the  western  Pacific.  Six  months  later  in  July 
(Figure  2.12a)  and  August  (Figure  2.12b)  the  reverse  is  true.  Towards  the  American 
coastline  the  trade  winds  have  a  more  north/south  velocity  component.  This  results  in 
a  more  sharply  defined  and  narrower  convergence  zone  between  the  two  trade  wind 
systems.  The  ITCZ  is  located  along  latitude  5°N  when  the  NE  trades  are  strongest 
during  the  Northern  Hemisphere  winter  (Figure  2.9)  and  is  positioned  at  latitudes  11- 
12 °N  when  the  SE  trade  winds  dominate  in  July  and  August  (Figure  2.12). 

To  sum  up,  the  high  pressure  center  over  the  eastern  Pacific  off  California  coupled 
with  the  intense  high  pressure  center  over  east  Asia  during  the  Northern  Hemisphere 
winter  cause  the  NE  trades  to  be  at  their  peak  and  dominate  the  SE  trade  wind  flow  in 
the  western  Pacific.  During  the  Northern  Hemisphere  summer,  the  high  center  west  of 
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California  moves  northwest,  away  from  the  equator  and  a  high  pressure  center  over  the 
south  Pacific  moves  north  toward  the  equator,  increasing  the  cross-equatorial  flow  of  the 
SE  trades  and  pushing  the  ITCZ  north  of  10°N  in  the  eastern  Pacific.  The  SE  trades 
cross  the  entire  length  of  the  equator  in  July  and  August.  In  the  western  Pacific,  the 
summer  heat  low  over  Asia  and  winter  cold  high  over  Australia  act  in  tandem  to  cause 
strong  southerly  flow  across  the  equator. 
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m.  RESULTS  AND  DISCUSSION 

A.      EDDY-RESOLVING  NATURE  OF  THE  GLOBAL  MODEL 

One  of  the  significant  features  of  the  Semtner  and  Chervin  general  circulation 
global  model  is  its  high  resolution  (1/2°  horizontal  grid  spacing  and  20  vertical  levels). 
This  enables  the  simulation  to  explicitly  resolve  mesoscale  eddies  in  the  tropics.  The 
Rossby  radius  of  deformation  is  the  critical  spatial  scale  numerical  models  need  to 
resolve  in  order  to  accurately  simulate  oceanic  boundary  phenomenon  such  as  fronts, 
eddies  and  boundary  currents.  The  Rossby  radius  varies  inversely  with  the  Coriolis 
parameter  /,  which  is  given  by 

/  =  1.47  x  HT'sinps'1, 
where  y?  is  latitude.  For  the  mid-latitude  oceans,  y?  =  45°,  the  Coriolis  parameter/  = 
1.0  x  10'V1.  At  10°  in  the  tropics  the  parameter/  =  0.25  x  10"4s_1.  This  factor  of 
four  difference  in  the  Coriolis  parameter  results  in  Rossby  radius  of  deformation  values 
typically  30-40  km  for  baroclinic  waves  in  mid-latitudes  and  120-160  km  in  the  tropics 
(Gill,  1982).  Thus  tropical  oceanic  phenomenon  should  be  accurately  modeled  by  the 
50  km  horizontal  resolution  of  the  SC-GCM.  The  model  also  uses  vertical  and  horizontal 
mixing  values  that  are  the  smallest  numerically  allowed  for  the  grid  spacing. 
Additionally  the  use  of  biharmonic  mixing  in  the  horizontal  allows  spontaneous  formation 
of  eddies  from  instabilities  without  the  undo  suppression  found  from  the  use  of  Laplacian 
mixing  in  other  model  formulations  (Semtner  and  Chervin,  1988). 
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The  near-surface  (37.5  meters)  instantaneous  current  vector  plots  (Figures  3. 1  and 
3.2)  of  15  May  and  14  November  in  the  model  show  a  dominant  train  of  anticyclonic 
eddies  across  the  central  and  eastern  Pacific  between  the  equator  and  5°N.  In  between 
these  eddies  and  to  the  south  can  be  seen  smaller  cyclonic  eddies.  The  anticyclonic 
eddies  appear  to  be  generated  from  the  horizontal  shear  zone  between  the  eastward 
flowing  North  Equatorial  Countercurrent  (NECC)  and  the  westward  flowing  South 
Equatorial  Current  (SEC).  The  modeled  anticyclonic  eddies  have  characteristic 
wavelengths  of  800-1000  km,  periods  of  20-30  days  and  propagate  westward  with  a 
phase  speed  of  40  km  day"1  (Tom  Murphree,  personal  communication).  They  are  most 
prominent  in  the  summer  and  fall  (Figure  3.2)  when  the  SE  trade  winds  are  at  their 
seasonal  maximum  and  the  modeled  NECC  and  SEC  reach  their  peak  in  response  to  this 
wind  forcing. 

The  observational  evidence  to  support  the  existence  of  these  eddies  is  abundant. 
Legeckis  (1977)  first  detected  westward  moving  long  waves  in  the  eastern  tropical  Pacific 
from  infrared  images  of  a  geostationary  satellite  in  1975.  The  satellite  images  displayed 
a  wavelike  temperature  front  between  the  SEC  and  NECC.  A  tongue  of  cooler  water 
projected  westward  by  the  SEC  results  from  the  advection  of  upwelled  water  along  the 
South  American  coast  as  well  as  upwelling  along  the  equator  from  the  divergence  of 
Ekman  transport.  The  NECC  advects  warmer  water  from  the  western  Pacific.  The 
satellite  analysis  of  these  long  waves  by  Legeckis  (1977)  revealed  wavelengths  of  800- 
1200  km,  periods  of  20-30  days  and  phase  speeds  of  40  km  day"1,  matching  the  SC- 
GCM  output.   The  waves  are  schematically  shown  in  Figure  3.3.    A  companion  study 
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of  the  instantaneous  thermal  field  in  the  SC-GCM  also  reveals  a  wavelike  temperature 
front  between  the  SEC  and  NECC  in  the  SC-GCM  (Tom  Murphree,  personal 
communication).  Interestingly,  Legeckis  noted  a  complete  absence  of  these  long  waves 
in  the  satellite  imagery  during  the  following  year,  1976,  and  related  it  to  anomalous 
warming  (4-6  °C)  in  the  eastern  Pacific  and  weaker  SE  trades  resulting  in  a  weakening 
of  the  SEC  and  NECC.  Further  analyses  of  satellite  imagery  from  1975-1981  by 
Legeckis  et  al.  (1983)  showed  these  long  waves  occurred  predominantly  during  the 
Northern  Hemisphere  summer  and  fall.  They  also  occurred  every  year  during  the  study 
except  1976  which  subsequent  research  revealed  was  an  El  Nino  year. 

Subsequent  studies  of  these  long  waves  with  satellite-tracked  drifter  buoys  (Hansen 
and  Paul,  1984)  and  comparisons  of  satellite  SST  with  ship  and  moored  buoy  data 
(Pullen,  et  al,  1987)  have  revealed  much  about  the  properties  of  these  waves.  Hansen 
and  Paul  (1984)  used  drifter  buoys  to  demonstrate  that  each  trough  in  the  temperature 
front  contains  an  anticyclonic  eddy  and  a  smaller  cyclonic  eddy  south  of  the  crest  of  the 
waves.  As  can  be  seen  in  Figure  3.2  of  the  SC-GCM  results,  the  model  accurately 
simulates  both  of  these  eddies. 

One  difference  noted  between  the  model  output  and  observations  is  that  the  model 
exhibits  the  anticyclonic  eddies  as  far  west  as  165  °E.  Evidence  from  observations 
(Wyrtlri,  1978)  indicate  that  the  furthest  westward  extent  yet  observed  in  these  waves  is 
159  °W.  Their  possible  existence  further  west  warrants  investigation  by  oceanographers. 

Barotropic  instability  from  the  shear  zone  between  the  SEC  and  NECC  is  believed 
to  play  a  major  role  in  the  formation  of  these  eddies  (Philander,  1978),  but  another 
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simulation  (Cox,  1980)  successfully  produced  these  waves  and  analysis  indicated  that 
baroclinic  instability  to  also  be  an  important  secondary  process.  Furthermore,  a  model 
study  by  Semtner  and  Holland  (1980)  produced  20-30  day  waves  without  a  NECC, 
indicating  that  the  role  of  barotropic  instability  in  their  formation  is  not  certain.  Further 
studies  of  the  results  in  the  SC-GCM  are  required. 

The  anticyclonic  mesoscale  eddies  are  easily  resolved  by  upwards  of  15-20  grid 
points  in  the  SC-GCM.  The  SC-GCM  additionally  resolves  the  smaller  cyclonic  eddy 
features  of  wavelengths  100-200  km.  A  poorer  resolution  model  would  not  be  able  to 
resolve  these  smaller  features. 

Mesoscale  eddies  have  a  critical  influence  on  the  mean  circulation.  They  transport 
momentum  and  heat  toward  parts  of  the  mean  circulation  and  away  from  other  parts. 
We  are  just  beginning  to  understand  the  geographic  distribution,  size  and  dynamics  of 
these  mesoscale  eddies.  In  parts  of  the  equatorial  Pacific,  mesoscale  eddies  are  so 
prevalent  and  strong,  as  in  Figure  3.2,  as  to  nearly  totally  obscure  the  mean  flow. 
Nevertheless  a  qualitative  description  of  the  seasonal  cycle  in  the  tropical  Pacific  current 
system  is  evident  in  the  vector  plots.  In  the  spring  (Figure  3.1)  the  Equatorial 
Undercurrent  (EUC)  attains  its  maximum  strength  and  can  be  seen  as  a  dominant 
eastward  flow  along  the  equator  at  the  37.5  meter  depth  in  the  model.  In  this  figure  it 
appears  that  the  EUC  is  gaining  strength  on  its  eastward  journey.  In  actuality  the  EUC 
and  its  high  velocity  core  shoals  towards  the  east;  i.e. ,  higher  velocity  isotachs  are  nearer 
the  surface  in  the  east  than  in  the  west.  Six  months  later  in  November  (Figure  3.2),  the 
EUC  has  weakened  and  none  of  its  isotachs  is  seen  at  the  37.5  meter  level  in  the  model 
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with  the  exception  of  eastward  velocity  on  the  equator  between  160°  and  170°E.  In  the 
late  summer  and  fall,  the  southeasterly  trades  attain  their  maximum  strength,  and  in 
response  the  SEC  reaches  maximum  velocity  and  transport.  The  NECC  similarly  is  at 
its  peak  magnitude  and  it  is  at  this  time  that  the  modeled  eddies  are  most  prominent. 

The  extremely  complicated  nature  of  oceanic  flow,  particularly  in  highly  variable 
regions  such  as  the  tropics  is  evident  in  these  plots  of  instantaneous  current  vectors. 
Shipboard  observations  have  difficulty  in  resolving  the  time  and  spatial  scales  of  these 
types  of  flows.  The  immense  expanse  of  water  to  be  surveyed  in  the  tropical  Pacific  and 
the  high  cost  of  ship  time  make  numerical  modeling  of  this  region  an  important  adjunct 
to  oceanographic  surveys  in  understanding  the  tropical  oceanic  circulation. 

Another  method  of  showing  the  instabilities  of  tropical  currents  is  with  time  series 
graphs  of  temperature.  Figures  3.4-3.7  begin  with  the  last  2.5  years  of  the  initial  model 
experiment  which  used  annual  mean  wind  forcing.  At  the  2.5  year  point  in  these  figures 
the  monthly  mean  wind  forcing  run  begins.  Average  temperature  from  20°S  to  20°N  at 
200  meters  is  plotted  at  various  longitudes  across  the  Pacific.  The  top  portion  of  each 
figure  is  a  3-day  running  mean  which  clearly  indicates  the  intrinsic  variability  of  the 
mean  equatorial  currents.  Eddy  activity  is  evident  in  both  the  annual  mean  and  seasonal 
mean  forcing  portions  of  the  graph  with  the  strongest  activity  in  the  central  Pacific  along 
140°W  (Figure  3.6).  The  beginning  of  the  seasonal  cycle  run  at  2.5  years  is  quite 
evident.  The  steady  rise  in  temperature  during  the  last  2.5  years  of  the  annual  run  is 
from  a  lack  of  intermediate-depth  convection  at  higher  latitudes.  This  warming  trend  is 
checked  with  the  introduction  of  forcing  by  monthly  varying  winds.  The  seasonal  cycle 
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run  quickly  reaches  a  quasi-equilibrium  state.  In  the  lower  row  of  plots,  high  frequency 
variability  is  filtered  out.  These  plots  are  180-day  and  360-day  running  mean 
temperature  values.  The  180-day  averages  clearly  show  a  semi-annual  signature  in  all 
the  time  series  plots  with  minimum  temperatures  during  the  Northern  Hemisphere 
summer  and  maximum  in  the  winter.  The  summer  temperatures  are  minimum  as  the 
upwelling  along  the  equator  as  well  as  between  the  NEC  and  NECC  are  at  their  annual 
maximum. 

B.      MEAN  ANNUAL  STRUCTURE  IN  THE  CENTRAL  PACIFIC 

The  seasonal  cycle  in  the  model  of  the  tropical  Pacific  is  analyzed  in  this  study 
from  archived  monthly  averaged  tapes  of  standard  oceanographic  field  variables  of  zonal 
velocity  (u),  temperature  (T)  and  salinity  (S).  The  analysis  of  horizontal  and  vertical 
sections  of  these  variables  is  simplified  by  using  these  tapes.  Although  the  model 
appears  to  resolve  quite  realistically  the  complex  structure  of  the  tropical  Pacific  in  the 
instantaneous  plots  as  demonstrated  in  the  previous  section,  a  validation  of  the  SC- 
GCM's  ability  to  simulate  the  seasonal  cycle  in  the  equatorial  Pacific  by  comparison  with 
observations  requires  an  analysis  of  averaged  tapes.  Because  of  the  large  annual  and 
interannual  variability  of  the  wind  field  in  the  Pacific,  a  comparison  of  instantaneous 
fields  from  a  climatologically  forced  model  with  observations  would  be  misleading.  The 
approach  of  this  study  is  to  use  monthly  averaged  tapes  and  make  comparisons  with  long- 
term  observational  studies.  The  forcing  over  the  tropical  Pacific  while  these  observations 
were  taken,  while  obviously  never  climatological,  was  at  least  typical  and  not  anomalous, 
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e.g.,  a  non  El  Nino  year.  The  complex  eddy  structure  that  is  resolved  by  the  model 
becomes  averaged  out  from  monthly  mean  tapes  and  the  mean  current  flow  can  more 
readily  be  evaluated.  Further  averaging  of  the  twelve  monthly  tapes  for  a  given  year 
produced  an  annual  mean  tape.  Although  the  eddies  are  "filtered  out"  from  the  monthly 
and  annual  mean  fields,  their  critical  role  and  large  influence  in  driving  the  mean 
circulation  has  been  completed. 

This  section  analyzes  the  annual  mean  tape  for  year  four  of  the  model  run.  As  was 
seen  in  the  time  series  of  temperature  (Figures  3.4-3.7),  the  model  reaches  a  quasi- 
equilibrium  state  quickly  after  the  introduction  of  seasonal  forcing.  (Year  four  of  the 
seasonal  run  is  year  6.5-7.5  in  these  figures.)  Because  of  the  lack  of  interannual 
variability  in  the  forcing  of  the  model,  one  model  year  looks  very  similar  to  any  other 
year.   The  choice  of  year  four  in  the  model  for  analysis  was  arbitrary. 

An  analysis  of  the  annual  mean  field  in  the  model  will  provide  the  backdrop  against 
which  to  understand  the  seasonal  cycle  of  the  Pacific  equatorial  current  system  and  allow 
direct  comparison  of  the  output  to  results  published  in  the  observational  literature.  The 
mean  annual  near-surface  (37.5  meters)  tropical  Pacific  currents  in  the  model  are  shown 
in  Figure  3.8.  Nearly  all  of  the  eddies  and  "wavy"  structure  previously  seen  in  the 
instantaneous  plots  are  no  longer  present.  The  currents  in  the  mean  are  zonal.  From 
150°W  to  100°W,  the  well-known  pattern  of  alternating  eastward  and  westward  currents 
in  the  equatorial  regions  can  be  plainly  seen.  The  eastward  flowing  EUC  at  the  equator 
splits  the  flow  of  the  SEC  into  a  strong  jet-like  branch  between  the  equator  and  5°N  and 
a  weaker  southern  branch  from  2°S  to  6°S. 
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The  mean  annual  meridional  sections  of  zonal  velocity,  temperature  and  salinity  in 
the  model  at  155  °W  longitude  are  compared  to  the  mean  annual  sections  from  the 
Hawaii-to-Tahiti  Shuttle  Experiment  conducted  between  1979-1980  in  the  central  Pacific. 
The  sections  from  this  experiment  were  analyzed  by  Wyrtki  and  Kilonsky  (1984).  The 
model  section  and  the  observations  show  many  striking  similarities.  Figures  3.9-3.11 
show  a  side  by  side  comparison  between  the  model  and  the  observations  for  zonal 
velocity,  temperature  and  salinity. 

Figure  3.9  and  Table  3.1  (page  24)  will  be  used  to  summarize  the  current 
comparisons  between  the  model  and  the  Hawaii-to-Tahiti  Shuttle  Experiment.  The 
currents  in  both  the  observed  and  the  modeled  sections  (Figure  3.9)  are  very  similar  both 
qualitatively  and  quantitatively.  Both  sections  show  a  North  Equatorial  Current  (NEC) 
with  nearly  the  same  latitudinal  extent  and  with  their  5  cm  s"1  isotach  extending  in  depth 
near  300  meters.  As  in  the  observations,  there  appears  a  tendency  for  the  NEC  in  the 
model  to  be  split  into  two  sections,  a  moderate  flow  of  more  than  20  cm  s"1  between 
10°N  and  15 °N  and  a  weaker  flow  of  5-10  cm  s"1  between  16  and  20°N  in  the  model 
compared  to  between  16  and  18°N  in  the  observations.  The  transports  of  the  westward 
NEC  is  calculated  at  25.2  Sv  in  the  model  and  23.3  Sv  for  the  Shuttle  experiment. 

The  North  Equatorial  Countercurrent  (NECC)  in  both  the  model  and  in  the 
observations  extends  from  4°N  to  9°N.  This  eastward  current  has  a  higher  velocity  core 
in  the  observations  than  in  the  model  (>40  cm  s"1  versus  27.7  cm  s"1).  The  5  cm  s"1 
isotach  extends  past  the  400  meter  depth  in  the  observations,  and  in  the  model  this 
isotach  stops  very  near  400  meters.  This  eastward  flowing  current  north  of  the  equator 
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is  actually  composed  of  two  separate  flows;  the  surface  NECC  and  the  deeper  North 
Subsurface  Countercurrent  (NSCC).  Both  the  model  and  observations  show  that  these 
two  eastward  currents  are  connected,  in  the  observations  they  are  also  connected  to  the 
EUC.  Wyrtki  and  Kilonsky  (1984)  used  a  velocity  minimum  at  170  meters  depth  to 
separate  these  two  flows  for  purposes  of  evaluating  the  mean  zonal  transport  of  the  two 
currents.  Such  a  division  in  the  model  output  is  not  evident,  but  for  purposes  of 
comparison,  the  model  level  of  160  meters  was  chosen  to  separate  the  NECC  and  NSCC. 
The  transport  of  the  surface  flowing  NECC  is  19.8  Sv  in  the  observations  and  13.3  Sv 
in  the  model.  The  difference  can  largely  be  attributed  to  the  higher  core  speed  of  the 
NECC  in  the  observations. 

The  North  Subsurface  Countercurrent  in  both  the  model  and  the  observations 
extends  from  4°  to  6°N.  It  also  shows  a  mean  observed  zonal  transport  of  8.9  Sv  and 
a  transport  of  6.7  Sv  in  the  model. 

The  Equatorial  Undercurrent  for  model  and  observations  straddles  the  equator 
from  2°S  to  2°N.  The  major  difference  is  the  larger  vertical  extent  of  this  thin  ribbon 
of  eastward  flow  in  the  model.  The  model  shows  the  currents  extending  to  over  400 
meters,  whereas  the  observations  show  a  reversal  of  the  flow  below  280  meters  in  the 
form  of  a  westward  flowing  Equatorial  Intermediate  Current  (EIC).  The  model  also 
shows  an  EIC  below  the  Undercurrent  but  it  appears  highly  variable  in  time  and  space. 
The  core  velocity  of  the  EUC  is  over  90  cm  s"1  at  a  depth  of  130  meters  in  the 
observations  and  111  cm  s"1  at  a  depth  of  160  meters  in  the  model.  Profile  current 
meter  measurements  of  the  EUC  taken  during  the  Shuttle  Experiment  gave  a  mean 
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transport  for  the  EUC  of  32.3  Sv  and  a  maximum  speed  of  102  cm  s"1.  In  contrast,  the 
model  EUC,  because  of  its  much  larger  vertical  extent,  had  a  transport  of  47  Sv.  Other 
researchers  (e.g.,  Tchernia,  1980)  give  a  mean  transport  for  the  EUC  of  40  Sv,  a  value 
between  the  model  results  and  the  observations.  Wyrtki  and  Kilonsky  (1984)  state  that 
the  EUC  was  unusually  weak  during  the  one-year  period  of  their  analysis. 

The  South  Equatorial  Current  is  divided  into  3  separate  sections  in  the  analysis  by 
Wyrtki  and  Kilonsky  (1984).  We  will  follow  a  similar  strategy  in  the  model  for  purposes 
of  comparison. 

Section  1  is  the  surface  westward  flow  of  the  SEC  from  4°N  to  the  equator.  This 
branch  is  much  more  restricted  in  the  observations  than  in  the  model  because  of 
"sandwiching"  by  the  connecting  eastward  flows  of  the  NECC  and  the  EUC.  The 
maximum  velocity  in  the  Shuttle  Experiment  for  the  NECC  is  nearly  90  cm  s"1  but  only 
57.5  cm  s-1  in  the  model.  The  model,  however,  shows  a  much  deeper  vertical  extent 
to  this  branch  of  the  SEC  because  of  the  separation  between  the  NECC  and  the  EUC. 
The  5cm  s" l  isotach  in  the  model  extends  beyond  the  435  meter  depth  of  the  figure.  The 
model  transport  for  this  section  of  the  SEC  is  19  Sv  compared  to  15.1  Sv  for  the 
observations. 

Section  2  of  the  SEC  in  the  model  and  the  Hawaii-to-Tahiti  Shuttle  Experiment 
extends  from  the  equator  to  9°S.  The  5  cm  s-1  isotach  for  both  is  similar  in  depth  of 
approximately  300  meters.  The  maximum  velocity  is  42.6  cm  s-1  in  the  model  and  just 
over  40  cm  s'1  in  the  observations.  Transport  for  this  branch  is  26.5  Sv  for  the  Shuttle 
experiment  versus  31.5  Sv  for  the  model.  In  both  the  model  and  the  observations,  this 
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section  of  the  SEC  has  by  far  the  highest  transport  relative  to  the  other  two  westward 
sections. 

Section  3  of  the  SEC  is  the  westward  flow  south  of  9°S.  Embedded  within  this 
flow  is  the  eastward  flowing  South  Equatorial  Countercurrent.  The  transport  of  this 
section  of  the  SEC  is  13.4  Sv  in  the  observations  and  8.3  Sv  in  the  model.  The  model 
shows  additional  westward  flow  in  this  branch  from  20 °S  to  17°S  of  an  additional  8.3 
Sv. 

The  total  transport  of  the  South  Equatorial  Current  from  its  3  sections  from  17°S 
to  4°N  is  55.0  Sv  for  the  Shuttle  Experiment  and  58.8  Sv  for  the  same  latitudinal  extent 
in  the  model.  The  additional  flow  in  the  model  south  of  17°S  brings  the  total  flow  for 
the  model's  SEC  to  67. 1  Sv.  Thus  the  SEC  has  the  greatest  magnitude  of  annual  mean 
volume  transport  of  all  the  equatorial  Pacific  currents.  This  can  be  attributed  to  the 
stronger  and  more  extensive  (latitudinally)  SE  trades  as  compared  to  the  NE  trades. 

The  South  Equatorial  Countercurrent  (SECC)  exhibits  two  branches  of  flow,  a 
weak  branch  from  8°  to  10°S  and  an  even  weaker  branch  at  12°S  in  the  observations  and 
16° S  in  the  model.  The  magnitude  of  the  velocities  of  the  flow  compare  favorably. 
Transport  of  only  0.6  Sv  was  calculated  during  the  Shuttle  Experiment  and  1.2  Sv  for 
the  model's  SECC.  Although  the  South  Equatorial  Countercurrent  is  slightly  stronger 
in  the  model  than  in  the  observations,  both  exhibit  the  similar  split  into  two  eastward 
branches. 

The  South  Subsurface  Countercurrent  (SSCC)  is  weak  in  the  model.  With  speeds 
of  less  than  5  cm  s"1  and  flow  from  4-9°S  at  250-430  meters,  the  eastward  zonal 
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Table  3.1   COMPARISON  OF  MODEL  TO  SHUTTLE 

EXPERIMENT 


MEAN  ANNUAL  CURRENTS 

MAXIMUM  ZONAL 
VELOCITIES  (cm  s"1) 

ZONAL  TRANSPORT  (Sv) 

CURRENT 

SHUTTLE 
EX. 

MODEL 

SHUTTLE 
EX. 

MODEL 

NEC 

20 

20.6 

23.3 

25.2 

NECC 

40 

27.7 

19.8 

13.3 

SEC  (eq.  to  4°N) 

90 

57.5 

15.1 

19.0 

SEC  (eq.  to  9°S) 

40 

42.6 

26.5 

31.5 

SEC  (9°  to  17°S) 

5 

9.0 

13.4 

8.3 

SEC  (Total) 

90 

57.5 

55.0 

58.8 

EUC 

102 

111.0 

32.3 

47.4 

NSCC 

10 

21.8 

8.9 

6.7 

SSCC 

5 

5.0 

4.3 

2.6 

SECC 

5 

9.0 

.6 

1.2 

transport  is  2.6  Sv.  The  observations  show  a  much  stronger  flow  with  a  transport  of  4.3 
Sv.  However,  both  the  model  and  the  observations  show  agreement  with  Tsuchiya 
(1975)  in  the  relative  strength  of  the  South  to  the  North  Subsurface  Countercurrent,  i.e., 
the  SSCC  is  half  as  intense  as  the  NSCC. 

Wyrtki  and  Kilonsky  (1984)  also  noted  an  unnamed  westward  flow  below  the 
NECC  between  6°  and  7°N  that  had  a  mean  transport  of  2.5  Sv.  The  model  resolves 
a  remarkably  similar  flow  in  the  same  area  with  a  westward  zonal  transport  of  2.7  Sv. 
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A  comparison  of  model  and  observed  temperature  structure  (Figure  3.10)  shows 
the  trough-ridge  system  that  is  a  well-observed  feature  of  the  equatorial  Pacific.  The 
modeled  thermocline  comes  very  close  to  the  surface  at  10  °N  due  to  the  up  welling 
caused  by  the  strong  Ekman  divergence  between  the  eastward  flowing  NECC  and  the 
westward  flowing  NEC.  The  thermocline  also  shows  the  characteristic  spreading  of  the 
isotherms  at  the  equator.  This  is  associated  with  the  EUC.  Upwelling  takes  place  above 
the  EUC  velocity  core  and  downwelling  below  the  core.  The  thermocline  in  the  model 
does  not  have  as  sharp  a  vertical  gradient  of  temperature  as  in  the  observations,  but  both 
observations  and  model  show  the  thermocline  beginning  to  deepen  and  spread  as  vertical 
temperature  gradients  decrease  towards  the  poles  north  of  10 °N  and  south  of  3°S.  The 
pools  of  warm  water  (27-28 °C)  at  10° N  and  10 °S  are  associated  with  the  advection  of 
warm  water  from  the  western  Pacific  by  the  eastward  flowing  NECC  and  SECC. 

Close  examination  of  the  observed  temperature  section  shows  that,  as  Wyrtki  and 
Kilonsky  (1984)  stated,  the  isotherms  below  the  lower  boundary  of  the  thermocline  are 
bent  downward  close  to  the  equator,  an  indicator  of  the  horizontal  temperature  gradients 
associated  with  the  North  and  South  Subsurface  Countercurrents  at  5°N  and  5°S.  The 
model  shows  a  similar  bending  of  the  10°  and  1 1  °C  isotherms  but  evidently  not  a  strong 
enough  horizontal  gradient  to  generate  a  strong  South  Subsurface  Countercurrent. 

The  salinity  distribution  (Figure  3.11)  shows  many  similarities  between 
observations  and  simulations.  The  salinity  increases  southward  across  the  equator  until 
a  high  salinity  cell  is  reached  in  the  subtropical  South  Pacific.  The  subsurface  salinity 
maximum  is  36. 1  psu  in  the  thermocline  and  is  very  nearly  identical  to  the  observations. 
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Both  sections  show  a  tongue  of  high  salinity  water  in  the  southern  tropics  penetrating  to 
the  equator  and  into  the  EUC.  There  is  another  salinity  maximum  in  the  thermocline 
north  of  the  equator.  It  penetrates  southward  to  about  15 °N  in  the  model  but  11°N  in 
the  observations.  Both  sections  exhibit  the  salinity  minimum  in  the  northern  tropics 
rising  up  from  a  depth  of  over  400  meters  at  20° N  to  a  depth  of  100  meters  at  10°N. 

C.      SEASONAL  CYCLE  IN  THE  EQUATORIAL  PACIFIC 

The  seasonal  variation  of  the  surface  winds  over  the  equatorial  Pacific  Ocean 
determine  the  location  of  the  Intertropical  Convergence  Zone  (ITCZ),  the  area  of  light 
winds  between  the  converging  trade  wind  systems  of  the  Northern  and  Southern 
Hemisphere.  When  the  equatorial  trade  winds  are  weak  during  the  Northern  Hemisphere 
spring,  the  ITCZ  is  near  the  equator.  The  surface  currents  are  weak  and  the  Equatorial 
Undercurrent  (EUC)  is  at  its  peak.  Figure  3. 12  shows  near-surface  (37.5  meters)  current 
vectors  in  the  east  central  Pacific  during  the  month  of  May.  Most  of  the  30-day  and 
shorter  period  waves  have  been  averaged  out.  The  EUC  dominates  the  flow  with  its  high 
velocity  jet  along  the  equator. 

Figure  3.13  depicts  average  currents  during  the  month  of  November.  The  ITCZ 
has  shifted  northward  to  near  10 °N  in  the  eastern  Pacific.  The  southeast  trades  are 
strong  and  the  surface  currents  respond  by  reaching  maximum  velocities.  The  SEC 
attains  speeds  in  excess  of  90  cm  s-1  and  a  latitudinal  extent  of  5°S-5°N.  The  NECC 
can  be  followed  unbroken  across  the  entire  tropical  Pacific,  achieving  velocities  of  over 
65  cm  s"1.    The  NEC  continues  to  remain  a  relatively  weak  flow  with  only  a  minor 
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seasonal  signature  but  does  appear  better  organized  than  in  the  spring.  Its  velocities 
approach  20-30  cm  s"1.  The  EUC  does  not  surface  during  November  and  is  significantly 
weaker  than  in  the  spring  with  a  maximum  velocity  of  118  cm  s"1  at  131°W. 

A  meridional  section  along  160°W  (Figure  3.14)  shows  zonal  velocity  from  20°S 
to  20  °N  for  the  months  of  May  and  November.  The  EUC  decreases  in  velocity  from 
121  cm  s"1  to  96.5  cm  s"1.  The  depth  of  its  maximum  velocity  core  changes  from  1 17.5 
meters  to  160.0  meters  as  it  weakens  and  submerges.  In  contrast  the  North  Equatorial 
Countercurrent  nearly  doubles  in  velocity  from  25  cm  s'1  to  49  cm  s"1.  Likewise,  the 
South  Equatorial  Current  increases  from  36.0  cm  s-1  in  the  spring  to  62. 1  cm  s-1  in  the 
fall.  The  North  Equatorial  Current  does  not  appear  to  have  a  strong  seasonal  fluctuation, 
both  the  spring  and  the  fall  velocities  are  in  the  neighborhood  of  20  cm  s~l.  The 
eastward  flowing  South  Equatorial  Countercurrent  shows  a  similar  pattern  to  its  northern 
counterpart.  It  is  a  weak  4  cm  s" 1  in  the  spring  and  increases  significantly  in  velocity 
in  the  fall  to  over  15  cm  s"1. 

The  temperature  structure  exhibits  a  strong  seasonal  signature  in  the  area  between 
the  North  Equatorial  Countercurrent  at  9°N  and  the  branch  of  the  South  Equatorial 
Current  north  of  the  equator.  In  the  spring  the  thermocline  has  a  weak  trough  near  3°N 
and  a  slight  ridge  at  10°N.  Latitudinal  temperature  gradients  are  relatively  weak.  The 
situation  is  reversed  in  the  fall.  The  meridional  temperature  sections  in  Figure  3.15 
show  this  remarkable  seasonal  change  to  the  isotherm  pattern.  The  well-observed  pattern 
of  isothermal  spreading  of  the  thermocline  from  the  influence  of  the  EUC  is  also  resolved 
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well  in  the  model.  The  isotherms  in  the  fall  are  pushed  upward  considerably  by  the 
upwelling  at  10°N,  creating  a  strong  horizontal  temperature  gradient  across  the  front. 

The  seasonal  variation  in  salinity  structure  of  the  model  at  160°W  is  quite  small 
(Figure  3. 16).  The  front  that  exists  in  the  fall  between  the  NECC  and  the  SEC  generates 
a  salinity  front  in  addition  to  the  temperature  front  previously  mentioned.  This  can  be 
seen  in  the  vertical  isohalines  between  5°N  and  10°N  which  separates  the  salty  water 
along  the  equator  from  the  slightly  fresher  water  to  the  north. 

Another  view  of  the  seasonal  changes  in  current  structure  along  the  equator  is 
shown  in  Figure  3.17.  This  represents  zonal  velocity  in  the  model  as  an  equatorial 
section  to  a  depth  of  310.0  meters.  Several  striking  features  stand  out.  The  top  figure 
of  flow  in  May  shows  that  the  EUC  reaches  a  maximum  strength  of  148  cm  s"1  at 
130°W  at  a  depth  of  1 17.5  meters.  There  is  eastward  flow  at  nearly  all  depths  along  the 
equator  as  far  west  as  140°E  where  the  current  is  very  deep  (222.5  meters),  but  it 
appears  that  the  EUC  is  a  poorly  organized  flow  until  it  approaches  the  international 
dateline.  During  the  spring  months  the  EUC  surfaces  and  there  is  eastward  flow  along 
the  surface  from  175  °W  to  near  the  South  American  continent.  In  November  the 
modeled  EUC  shows  a  maximum  velocity  of  118  cm  s"1,  again  along  130°W  longitude 
but  at  a  deeper  160  meters  depth.  At  this  time  of  the  year,  the  flow  at  the  surface  along 
the  equator  is  westward  under  the  influence  of  the  SE  trades.  Another  feature  of  the 
model's  EUC  is  the  pronounced  upward  slope  toward  the  east  of  the  velocity  core.  This 
is  characteristic  of  the  observations  of  the  EUC  as  it  straddles  the  upward  sloping 
thermocline  on  its  eastward  journey  across  the  Pacific. 
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Figure  3.18  is  a  time-longitude  plot  at  the  first  vertical  level  of  the  model  (12.5 
meters)  that  shows  zonal  velocity  along  the  equator  as  a  function  of  time.  The  surfacing 
of  the  EUC  is  first  seen  in  the  spring  commencing  in  the  eastern  Pacific. 

The  seasonal  fluctuations  in  zonal  volume  transport  in  Sverdrups  for  the  NECC, 
EUC  and  SEC  at  various  longitudes  are  shown  in  Figure  3. 19.  The  SEC  and  the  NECC 
are  in  phase  with  each  other  and  out  of  phase  with  the  EUC  as  previously  mentioned. 
Figure  3.20  shows  the  variability  in  the  transport  of  the  NEC.  No  coherent  seasonal 
cycle  is  evident  with  this  westward  flow. 

D.      WESTERN  EQUATORIAL  PACIFIC 

The  currents  in  the  western  Pacific  are  of  an  extremely  complex  nature  due  to 
monsoon  cycle  forcing  which  causes  a  reversal  of  the  wind  stress,  and  due  to  the 
complex  topography  represented  by  the  hundreds  of  islands  of  the  Philippines  and  of  the 
Indonesian  Archipelago.  Many  of  the  equatorial  Pacific  currents  terminate  or  originate 
in  this  region  and  the  treatment  of  the  islands  in  numerical  models  determines  to  a  great 
extent  the  accuracy  of  the  reproduced  currents.  Figure  3.21  represents  the  observed 
current  structure  during  February  and  August  in  the  western  Pacific.  The  NEC 
bifurcates  along  the  Philippine  island  of  Mindanao  with  a  portion  of  its  flow  heading 
northwest  to  form  the  beginnings  of  the  Kurishio  and  the  remainder  of  its  transport 
deflecting  to  form  the  start  of  the  NECC.  The  SEC  joins  the  surfaced  New  Guinea 
Undercurrent  on  the  north  side  of  that  island  and  is  turned  by  the  island  of  Halmahera 
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to  join  in  forming  the  beginning  of  the  NEC.    The  importance  of  Halmahera  in  the 
current  structure  is  seen  in  the  historical  ship  drift  observations  in  Figure  3.22. 

The  SC-GCM  in  its  global  formulation  was  forced  to  make  compromises  with  the 
global  topography  in  order  to  make  the  global  calculation  tractable.  Many  of  these 
topographic  compromises  were  made  in  the  western  Pacific  because  of  the  large  number 
of  islands  there.  For  example,  the  islands  of  Mindanao  and  Halmahera  were  shoaled. 
This  results  in  model  currents  for  the  western  Pacific  seen  in  the  vector  plot  of  Figure 
3.23.  Since  Mindanao  is  shoaled,  the  NEC  flows  to  the  coast  of  Borneo  (5°N,  120°E) 
and  is  shunted  toward  the  south  to  form  an  intense  unstable  current  along  the  coast  in  the 
Makassar  strait.  The  large  anticyclonic  eddy  just  north  of  the  equator  between  120°E 
and  130°E  becomes  a  semi-permanent  feature  in  the  model  calculation  and  represents  the 
origin  of  the  eastward  flowing  NECC.  The  fact  that  the  central  Pacific  currents  in  the 
model  examined  earlier  compared  so  well  to  observations  indicates  that  the  treatment  of 
the  islands  in  the  west  and  the  boundary  they  represent  are  quickly  "forgotten"  by  the 
currents  after  having  reached  the  central  Pacific  region. 
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IV.  CONCLUSIONS  AND  RECOMMENDATIONS 

The  Semtner  and  Chervin  global  ocean  model  reproduces  many  of  the  known 
features  of  the  seasonal  cycle  of  the  tropical  Pacific  Ocean.  Current  structure,  velocity 
and  seasonal  variability  show  remarkably  realistic  features.  However,  many  of  the 
details  of  these  plots  have  no  corresponding  observations  with  which  to  make  direct 
comparisons.  The  observations  that  are  available  often  cannot  be  meaningfully  compared 
to  model  output  because  of  the  difference  in  wind  forcing  at  the  time  observations  were 
taken  from  the  climatological  wind  forcing  in  the  model. 

Temperature  and  zonal  velocity  patterns  fit  the  descriptions  of  these  fields  in  the 
observational  literature.  One  notable  exception  is  that  the  observed  mixed  layer  depths, 
particularly  in  the  western  Pacific,  are  deeper  than  in  the  model.  The  model  tends  to 
have  a  diffuse  thermocline.  Philander  et  al.  (1987)  noted  very  similar  results  from  his 
model  studies  and  believed  that  the  poor  mixing  might  be  attributed  to  the  lack  of  high- 
frequency  wind  fluctuations  in  the  climatological  wind  forcing  fields. 

A  detailed  analysis  of  the  mass  and  heat  budgets  from  the  SC-GCM  needs  to  be 
made  to  ascertain  the  model's  performance  in  these  areas.  Similar  analyses  and 
comparisons  with  the  results  of  the  model  in  other  tropical  regions  and  with  other 
regional  models  should  also  be  made.  Analysis  of  the  deep  equatorial  flow  would  be  of 
interest  to  determine  if  the  off-equatorial  and  deep  equatorial  jets  that  have  been  observed 
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are  resolved.  These  deep  equatorial  jets  have  also  appeared  in  model  studies  by 
McCreary  (1987). 

The  complexity  of  oceanic  flow  in  the  western  Pacific  requires  a  more  sophisticated 
treatment  of  islands  in  the  model  in  order  for  this  area  to  show  more  realism.  Kill  worth 
etal.  (1991)  have  developed  a  free-surface  effect  that  should  make  calculations  involving 
islands  tractable. 

Improvement  in  the  horizontal  and  vertical  resolution  of  the  model,  while  desirable, 
would  not  likely  improve  the  model  results  as  much  as  improving  the  climatological 
forcing  of  wind  and  heat.  A  test  of  the  model  using  observed  wind  forcing,  particularly 
incorporating  interannual  variability,  would  be  of  great  interest. 


32 


APPENDIX  (FIGURES) 


Figure  1.1   A  schematic  diagram  of  upper  tropical  Pacific  Ocean  circulation  (from 
Philander  et  al. ,  1987) 
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Figure  2.1  The  global  ocean  model  domain  for  the  Pacific  Ocean. 
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Figure  2.3  Hellerman  and  Rosenstein  (1983)  mean  wind  stress  in  western  and  central 
Pacific  for:  (a)  January  and  (b)  February.   Standard  vector  length  =  1  dyn  cm"2. 
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Figure  2.4  Hellerman  and  Rosenstein  (1983)  mean  wind  stress  in  western  and  central 
Pacific  for:  (a)  March  and  (b)  April.   Standard  vector  length  =  1  dyn  cm-2. 
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Figure  2.5  Hellerman  and  Rosenstein  (1983)  mean  wind  stress  in  western  and  central 
Pacific  for:  (a)  May  and  (b)  June.   Standard  vector  length  =  1  dyn  cm"2. 
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Figure  2.6  Hellerman  and  Rosenstein  (1983)  mean  wind  stress  in  western  and  central 
Pacific  for:  (a)  July  and  (b)  August.   Standard  vector  length  =  1  dyn  cm"2. 
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Figure  2.7  Hellerman  and  Rosenstein  (1983)  mean  wind  stress  in  western  and  central 
Pacific  for:  (a)  September  and  (b)  October.   Standard  vector  length =1  dyn  cm"2. 
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Figure  2.8  Hellerman  and  Rosenstein  (1983)  mean  wind  stress  in  western  and  central 
Pacific  for:  (a)  November  and  (b)  December.   Standard  vector  length  =  1  dyn  cm"2. 
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Figure  2.9  Hellerman  and  Rosenstein  (1983)  mean  wind  stress  in  eastern  Pacific  for: 
(a)  January  and  (b)  February.   Standard  vector  length  =  1  dyn  cm"2. 
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Figure  2.10   Hellerman  and  Rosenstein  (1983)  mean  wind  stress  in  eastern  Pacific 
for:  (a)  March  and  (b)  April.   Standard  vector  length =1  dyn  cm"2. 
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Figure  2.11    Hellerman  and  Rosenstein  (1983)  mean  wind  stress  in  eastern  Pacific 
for:  (a)  May  and  (b)  June.   Standard  vector  length  =  1  dyn  cm"2. 
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Figure  2.12   Hellerman  and  Rosenstein  (1983)  mean  wind  stress  in  eastern  Pacific 
for:  (a)  July  and  (b)  August.   Standard  vector  length  =  1  dyn  cm"2. 
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Figure  2.13   Hellerman  and  Rosenstein  (1983)  mean  wind  stress  in  eastern  Pacific 
for:  (a)  September  and  (b)  October.   Standard  vector  length =1  dyn  cm"2. 
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Figure  2.14   Hellerman  and  Rosenstein  (1983)  mean  wind  stress  in  eastern  Pacific 
for:  (a)  November  and  (b)  December.   Standard  vector  length =1  dyn  cm"2. 
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Figure  3.1    Instantaneous  current  vectors  at  37.5  meters  for  15  May,  year 
seasonal  model  run.   a)  Pacific  basin  b)  Central  eastern  Pacific. 
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Figure  3.2  Instantaneous  current  vectors  at  37.5  meters  for  14  November,  year  4  of 
seasonal  model  run.   a)  Pacific  basin  b)  Central  eastern  Pacific. 
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Figure  3.3  Schematic  representation  of  wave-like  sea  surface  temperature  front  in  the 
eastern  equatorial  Pacific.  The  anticyclonic,  eddy-like  patterns  were  inferred  from 
drifting  buoy  tracks,  (from  Legeckis  et. al.,  1983). 
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Figure  3.4  Time  series  of  integrated  temperature  along  200  m  depth  from  20°N-20°S 
at  160°E:  (a)  3-day  running  mean,  (b)  180-day  running  mean,  and  (c)  360-day 
running  mean. 
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Figure  3.5  Time  series  of  integrated  temperature  along  200  m  depth  from  20°N-20°S 
at  170° W:  (a)  3-day  running  mean,  (b)  180-day  running  mean,  and  (c)  360-day 
running  mean. 


52 


M0ufsr  200M   20N-J0S  TEMP 


140UEST  200M   20N-20S  TEMP 


_i 1 i L 1 I i_ 


a  i  2  :  *  5  t> 


M0UEST     200M       20N-20S     IEMP 


i ' r 


.-A. 


V 


_L    —I J I L L. 


'  1  3 


■»  ia 


Figure  3.6  Time  series  of  integrated  temperature  along  200  m  depth  from  20°N-20°S 
at  140° W:  (a)  3-day  running  mean,  (b)  180-day  running  mean,  and  (c)  360-day 
running  mean. 
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Figure  3.7  Time  series  of  integrated  temperature  along  200  m  depth  from  20°N-20°S 
at  110°W:  (a)  3-day  running  mean,  (b)  180-day  running  mean,  and  (c)  360-day 
running  mean. 
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Figure  3.8  Mean  annual  near  surface  currents  (37.5  meters)  in  central  eastern  Pacific 
for  year  4. 
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Figure  3.9  Meridional  section  of  mean  annual  zonal  velocity  during  a)  Hawaii-to- 
Tahiti  Shuttle  Experiment  (from  Wyrtki  and  Kilonsky,  1984)  and  b)  year  4  of  the 
model.   Contour  interval  is  5  cm  s"1. 
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Figure  3.10  Meridional  section  of  mean  annual  temperature  distribution  during  a) 
Shuttle  Experiment  (from  Wyrtki  and  Kilonsky,  1984)  and  b)  year  4  of  the  model. 
Contour  interval  is  1°C  and  depth  limit  is  435  meters  in  the  model  figure. 
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Figure  3.11  Meridional  section  of  mean  annual  salinity  distribution  during  a)  Shuttle 
Experiment  (from  Wyrtki  and  Kilonsky,  1984)  and  b)  year  4  of  the  model.  Contour 
interval  is  .  1  psu  and  depth  limit  is  435  meters  in  the  model  figure. 
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Figure  3.12  Mean  monthly  current  vectors  for  May  of  Year  4  in  central  Pacific 
portion  of  the  model .  Surfacing  of  Equatorial  Undercurrent  is  apparent  along  equator. 
Depth  depicted  is  37.5  meters. 
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Figure  3.13  Mean  monthly  current  vectors  for  November  of  Year  4  in  central 
Pacific.  Equatorial  Undercurrent  has  submerged  and  North  Equatorial  Countercurrent 
and  South  Equatorial  Current  are  now  dominant.   Depth  depicted  is  37.5  meters. 
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Figure  3.14  Meridional  section  along  160°W  of  mean  monthly  zonal  velocity  during 
a)  May  of  year  4  and  b)  November  of  year  4.   Contour  interval  is  5  cm  s"1. 
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Figure  3.15  Meridional  section  along  160°W  of  mean  monthly  temperature 
distribution  during  a)  May  of  year  4  and  b)  November  of  year  4.  Contour  interval 
is  1°C. 
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Figure  3.16  Meridional  section  along  160°W  of  mean  monthly  salinity  distribution 
during  a)  May  and  b)  November  of  year  4.   Contour  interval  is  .2  psu. 
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Figure  3.17    Equatorial  section  of  zonal  velocity  during  a)  May  of  year  4  and  b) 
November  of  year  4.   Contour  interval  is  10  cm  s"1. 
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Figure  3.18   Longitude-time  plot  of  zonal  velocity  along  the  equator  during  Year  4 
at  12.5  meters.   Contour  interval  is  10  cm  s"1. 
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Figure  3.19  Zonal  volume  transport  of  mean  monthly  equatorial  currents  along 
various  meridians,  a)  North  Equatorial  Countercurrent  b)  Equatorial  Undercurrent  c) 
South  Equatorial  Current. 
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Figure  3.20  Zonal  transport  of  North  Equatorial  Current  in  Sverdrups  (106  m3  s"1). 
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Figure  3.21  Current  pattern  in  the  western  Pacific  during  a)  February  and  b)  August, 
(from  Tchernia,  1980). 
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Figure  3.22   Historical  monthly-mean  ship  drift  observations  in  the  western  Pacific 
from  May  through  November,  (from  Hacker,  et.al.,  1989) 
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Figure  3.23  Instantaneous  current  vectors  in  western  Pacific  for  14  November,  year 
4  of  model. 
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